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Photoionization of hydrogen atom by coherent intense high-frequency short laser 
pulses: Direct propagation of electron wave packets on enormous spatial grids 
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The time-dependent Schrodinger equation for the hydrogen atom and its interaction with coherent 
intense high-frequency short laser pulses is solved numerically exactly by employing the code imple- 
mented for the multi-configurational time-dependent Hartree-Fock (MCTDHF) method. Thereby, 
the wavefunction is followed in space and time for times longer than the pulse duration. Results 
are explicitly shown for 3 and 10 fs pulses. Particular attention is paid to identifying the effect 
of dynamic interference of photoelectrons emitted with the same kinetic energy at different times 
during the rising and falling sides of the pulse predicted in [Ph. V. Demekhin and L.S. Cederbaum, 
Phys. Rev. Lett. 108, 253001 (2012)]. In order to be able to see the dynamic interference pattern 
in the computed electron spectra, the photoelectron wave packet has to be propagated over long 
distances. Clearly, complex absorption potentials often employed to compute spectra of emitted 
particles cannot be used to detect dynamic interference. For the considered high-frequency pulses 
of 3 and 10 fs durations, this requires enormously large spatial grids. The presently computed 
photoionization and above-threshold ionization spectra are found to exhibit pronounced dynamic 
interference patterns. Where available, the pattern is in a very good agreement with previously 
published results. 
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PACS numbers: 33.20.Xx, 41.60.Cr, 82.50.Kx 



I. INTRODUCTION 

The presently available light sources, like attosecond 
lasers [11] , high-order harmonic generation sources 0, Q ^ 
or free electron lasers 0, Qi allow one to study the 
interaction of matter with super-intense coherent high- 
frequency ultra-short laser pulses. Many new phenom- 
ena, which are not available or difficult to observe with 
optical lasers operating in the nano- and picosecond 
regimes, arise in experiments with such pulses p-lOf. One 
of the main advantages here is that the high carrier fre- 
quencies allow one to directly access a few well-separated 
highly-excited electronic states of a system. This is usu- 
ally not possible with optical pulses triggering simulta- 
neously a dense spectrum of electronic states. Another 
advantage is that short pulses allow one to study the 
initiated dynamics on the relaxation timescale of highly- 
excited states (typically femto- or even attoseconds). The 
theoretical description of the processes triggered by re- 
alistic pulses requires the propagation of electron wave 
packets created by the pulse in real time and real. 

One of the phenomena arising due to the interaction of 
matter with coherent intense high-frequency short pulses 
is known as dynamic interference [lO[. There, the pulse 
ionizes a system by the absorption of a single photon 
and; at the same time, induces a time-dependent shift 
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of the 'dressed' ground state of an atom relative to the 
continuum due to an AC Stark effect in the electronic 
continuum. The energy shift adiabatically follows the 
pulse intensity envelope g^ (t) [Il[ . Because of this time- 
dependent energy shift, the photoelectrons emitted when 
the pulse rises have the same kinetic energy as those emit- 
ted when the pulse decreases. These two electron wave 
packets emitted at different times superimpose and inter- 
fere, and the resulting photoelectron spectrum exhibits 
a pronounced multipeak pattern [lO[. Similar dynamic 
interference of electrons was also found in the above- 
threshold ionization spectra of model anions [IJ, [l^ , in 
the sequential multiphoton ionization of atoms exposed 
to optical [IJ] and high-frequency [i5| pulses, as well as 
in the resonant Auger effect of atoms in free electron laser 
pulses |15| . 

In order to solve the time-dependent Schrodinger equa- 
tion for the hydrogen atom interacting with a pulse we 
have used in our recent work [lO| a previously developed 
theoretical approach |16l - [2]| . In this approach the total 
wave function is expanded in terms of the full set of the 
field-free stationary states of the system. This leads to 
equations for the corresponding population amplitudes 
which were then propagated in a large but restricted rel- 
evant subset of these states. In the present work, we 
solve the time-dependent Schrodinger equation numer- 
ically exactly for the same problem by directly propa- 
gating the single-electron wave packet during and briefly 
after the pulse and extract from this wave packet the full 
spectrum of the emitted electrons. This is a challenging 
problem by itself (see below). To allow for the dynamic 



interference to take place, one needs to propagate the full 
photoelectron wave packet in space and time without us- 
ing complex absorption potentials at the boundaries. For 
high-frequency pulses of even a few femtosecond duration 
the propagation requires spatial grids of ^ W^ Bohr. We 
have found that cutting off even a small tail of the wave 
packet falsifies the results seemingly. 

To the best of our knowledge, no such explicit time- 
dependent calculations on the electron wave packet prop- 
agation are reported in the literature for smooth realistic 
pulses, even not for a single-electron system. There is 
a series of works |22h24| (see also references therein for 
earlier works by the same authors), reporting calcula- 
tions performed for a smoothed Coulomb potential and 
short trapezoid-shaped pulse envelopes using restricted 
spatial grids and a complex absorption potential at the 
boundary. The pulses used consist of a rapidly rising 
edge, a long plateau, and a rapidly falling edge "to let 
the major fraction of photoelectrons to be formed at the 
pulse plateau". We note that such a pulse cannot give 
rise to dynamic interference by definition even if a com- 
plex absorbing potential were not used. The calculations 
performed for w = 30 eV [2j, [23| showed that the pho- 
toelectron peak position is shifting to higher energies as 
the intensity of the field grows. This finding already in- 
dicates that the binding energy of the system decreases 
with the field strengths |23| . which to our opinion is due 
to the AC Stark effect in the electronic continuum. 

In order to solve the technically challenging prob- 
lem at hand one needs sophisticated techniques. The 
multi-configuration time-dependent Hartree (MCTDH) 
method [25|, [2a| is known to be an optimal approach for 
wave packet propagation in many degrees of freedom and 
this approach requires efficient propagation algorithms 
for the underlying single-particle time-dependent func- 
tions [231 ■ Being formulated for bosonic particles the 
method is known as MCTDHB |2g, and for fermions as 
MCTDHF [li-lH. In the present work we utilize a par- 
ticular code of MCTDHF [SjI which is based on the gen- 
eral formulation of the problem given in |3J] and has been 
slightly modified to efficiently propagate single-particle 
wave packets in the presence of strong pulses. 



II. THEORY 

The total Hamiltonian for the hydrogen atom inter- 
acting with the linearly polarized laser field is given by 
(atomic units are used throughout) 
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where the electric field 



£{t) ~ 8og{t)cosujt. 



(1) 



(2) 



is polarized along the z-axis. Here, Eq is the field ampli- 
tude and uj is the carrier frequency of the pulse with a 



time-envelope g{t). The approach we use here to solve 
the time-dependent Schrodinger equation for the Hamil- 
tonian ([T|) is described in detail in [SjI and thus only a 
few points of essential relevance to the present work are 
outlined below. 

First, we make use of the axial symmetry along the z 
axis of the problem and employ a partial wave expansion 
of the single-electron wave function ^(r,t) in terms of 
spherical harmonics Ygo 
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By substituting the wave function ^ into the time- 
dependent Schrodinger equation for the Hamiltonian ([1} , 
one straightforwardly obtains the following system of 
coupled equations for the radial harmonics Pe{r,t) 
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The radial coordinate in the above system of coupled 
one-dimensional equations ^ can further be represented 
by a discrete variable representation (DVR) basis set 



P,(r,t)=^6,,,(t)x.(r) 



(5) 



As usual, employing DVRs has the advantage that spa- 
tially local operators posses a diagonal representation 
where a function acquires the value on the respective 
DVR grid point [33 . In order to avoid full matrices of 
the kinetic energy operator, we use here the finite element 
discrete variable representation (FEDVR) introduced in 
[36| . We thus divide the radial coordinate space in a cho- 
sen number of finite elements. In each finite element, the 
basis functions Xi('') are represented by the normalized 
Legendre interpolating polynomials 
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constructed over a Gauss-Lobatto grid {r^} with weights 
{i(Ji}. Note that each basis function vanishes at all grid 
points except one. Correspondingly, one arrives at a 
banded structure of the kinetic energy matrix, which 
makes numerical solution of Eqs. ([4]) faster, and, at the 
same time, we can use sufficiently flexible basis sets ad- 
justable to different physical problems. 

As a consequence of the DVR representation (O, the 
system of equations (j4|) can be straightforwardly trans- 
formed to a system of equations for the evolution of the 
time-dependent expansion coefficients bi^i(t) with known 
one-particle integrals. This system was propagated using 



the short-iterative Lanczos method employing algorithm 
[37| to approximate the exponential time-evolution oper- 
ator. The initial ground state of the H atom in the ab- 
sence of the pulse (although analytically known) is consis- 
tently obtained via imaginary time propagation starting 
from a guess function, which requires a negligible time as 
compared to the real time propagation. 

The final momentum distribution of the emitted pho- 
toelectrons was computed by Fourier transforming the 
electron wave packet ^(r) = ^(r,i = oo) after the pulse 
has expired 
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In order to exclude contributions to the final wave packet 
^(r) from electrons remaining bound to the nucleus (i.e., 
from the ground and Rydberg states), the inner-part re- 
gion of the radial spatial variable r was excluded from 
the transformation ([7]). Finally, the kinetic energy spec- 
trum (j{£) of the emitted photoelectrons is obtained by 
angular-averaging the momentum density distribution as 



a{e) ^ k J \^{k)\^dn. 



(8) 



and using that k = \/2e. 

As in our previous study [lO[, the present calcula- 
tions were performed for a Gaussian-shaped pulses with 



a time-envelope g{t) 



-t'/r' 



and central frequency 



u! — 53.6057 eV, which is 40 eV above the ionization 
threshold of H. The electron wave packets were propa- 
gated in the time interval of [— 3t, -t-Sr], for which the 
field amplitude at the interval boundaries is almost four 
orders of magnitude weaker than at the pulse maximum 
(t = 0). Calculations were performed for two pulse dura- 
tions of T = 3 fs and r = 10 fs. The shorter pulses were, 
but the longer pulses were not studied in our previous 
work ^. 

As will become evident below, the computed electron 
spectra exhibit a clear sequence of above threshold ion- 
ization (ATI) peaks located around the multiphoton ab- 
sorption energies £g = n ■ oj — IP. For the largest 
field intensity considered here, the third ATI peak lo- 
cated at £q = 4w — IP is about five orders of magnitude 
weaker than the main photoionization peak located at 
Eq = w — IP (see discussion around Fig. [2]). Because of 
this fact, we have decided to choose the parameters of the 
present calculations such that the photoionization peak 
and the three subsequent ATI peaks in the final elec- 
tron energy spectra are described accurately. For this 
purpose, we had to include the £ < A harmonics in the 
partial wave expansion ([3]). 

For the photon energy oj used, the third ATI peak 
is formed by electrons with momentum k ~ 3.85 a.u. 
Consequently, we have chosen the size of the radial grid 
such that photoelectrons with fc < 4 a.u. do not hit 
the outward grid boundary during the whole propaga- 
tion time. During the r = 3 fs pulse, the photoelec- 
trons with momentum k ^ A a.u. may move off the 



nucleus by Umax = Brfc ^ 3000 a.u., which was chosen 
as the radial grid size. The latter was represented by 
1000 equidistant finite elements of size 3 a.u. Each fi- 
nite element was covered by 15 Gauss-Labbato points. 
For the r = 10 fs pulse, the radial grid parameters were: 
Rmax = QTk ^ 10000 a.u. represented by 3125 finite el- 
ements with the size of 3.2 a.u., and each finite element 
is covered by 16 Gauss-Labbato points. The convergence 
of the solution with respect to the spatial grid represen- 
tation and integration time step has been ensured. 



III. RESULTS AND DISCUSSION 

Figure [T] illustrates the time evolution of the total elec- 
tron wave packet computed for the hydrogen atom ex- 
posed to short pulses of carrier frequency u ~ 53.6057 cV. 
The upper panel shows the results for the r = 3 f s pulse 
and peak intensity of /q = 7 x 10^^ W/cm^. Before the 
pulse arrives, the electron is in the Is ground state of H 
which is very close to the nucleus (note the scale of the 
figure). At very early times, when the pulse arrives, the 
part of the wave packet which will contribute to the pho- 
toclectron spectrum is created around the nucleus and 
starts to propagate outwards. Before the maximum of 
the pulse has arrived (time fs corresponds to the pulse 
maximum, thick blue curve), the wave packet being con- 
tinuously pumped by the pulse does not posses a maxi- 
mum. A maximum in the radial electron density distri- 
bution starts to develop only after the pulse maximum 
has arrived. For the wave packet computed at t = 3 fs 
the maximum of the wave packet is around r — 275 a.u. 

As the pulse expires, the maximum becomes more and 
more pronounced and it moves further outwards. At t = 
6 fs this maximum is around r = 525 a.u. In the final 
wave packet computed at i = 3r = 9 fs, the maximum is 
located at about r = 750 a.u. This maximum is mainly 
responsible for the photoionization peak at ej = w — IP 
in the final energy spectrum computed via Eqs. ([7]) and 
([5|. The weak humps seen in the final wave packet at 
about r = 1150 a.u. and 1500 a.u. are mainly formed 
by the photoelectrons contributing to the first and the 
second ATI peaks in the final energy spectrum. The lower 
panel of Fig. [T] depicts the time evolution of the electron 
wave packet for the r = 10 fs pulse and peak intensity of 
/o = 1-5 X 10^^ W/cm^. The results are rather analogous 
to those described above for the r = 3 fs pulse except 
that the time and the spatial scales arc different. The 
detached electron now reaches much longer distances and 
the calculation has become even much more cumbersome. 

The energy distribution of the electrons emitted dur- 
ing the illumination of the hydrogen atom by coherent 
T = 3 fs laser pulses of different peak intensities are de- 
picted in Fig. [21 For transparency of the figure, the en- 
ergy range is restricted to show the photoionization peak 
and the two subsequent ATI peaks (note the logarithmic 
scale on the vertical axis). One can see that, as the field 
strength increases from the bottom to the top of the fig- 
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FIG. 1: (Color online) Time-evolution of the radial electron 
density (total electron wave packet, |^(r, f)p) computed for 
the hydrogen atom exposed to Gaussian-shaped short pulses 
of carrier frequency of oj = 53.6057 eV. Upper panel: The 
duration of the pulse is r = 3 fs, and the peak intensity 
7o = 7 X 10^*^ W/cm^. The wave packet shown for f = -9 fs 
(straight orange line near the coordinates origin) is the ground 
state wave packet. Note that the final wave packet computed 
at t = 9 fs spreads out to about r ~ 3000 a.u. Lower panel: 
The duration of the pulse is r = 10 fs, and the peak intensity 
lo = 1.5 X 10^® W/cm^. The wave packet shown for t = 
—30 fs is the ground state wave packet. Note that the final 
wave packet computed at f = 30 fs spreads out to about 
r ~ 10000 a.u. 
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FIG. 2: Computed energy distributions of the emitted photo- 
electrons after exposure of the H atom to coherent laser pulses 
of r = 3 fs duration, u = 53.6057 eV carrier frequency, and 
different peak intensities (indicated in each panel). Shown are 
the photoionization peak and the two subsequent ATI peaks 
produced by multiphoton absorption from the ground state 
(the next ATI peaks are weak and are not shown). The field- 
free multiphoton absorption energy positions Sg = n- uj — IP 
are indicated by vertical gray lines. The insets in the up- 
permost panel illustrate on an enlarged scale the similarity 
between the dynamic interference patterns in the photoion- 
ization and ATI spectra (see also text for details). 



ure, the photoionization peak and both ATI peaks are 
systematically shifted to higher electron energies. For 
reference, the field-free multiphoton absorption energy 
positions Eg = n • a; — IP are indicated by vertical lines. 

All peaks in the electron energy distribution exhibit 
pronounced multiple-peak structures which are due to 
the dynamic interference [T^. This fact is demonstrated 
in the insets in the uppermost panel for the largest inten- 
sity considered here. The left inset shows the photoion- 
ization peak on an enlarged scale, whereas the middle and 
right insets enlarge the first and the second ATI peaks, 
respectively. As can be seen from these insets, the dy- 
namic interference patterns are qualitatively similar for 
the photoionization peak and for the two subsequent ATI 
peaks, each possessing three well-resolved oscillations of 
the intensity (the vertical axis in all insets is in logarith- 
mic scale) . The maximum of the intensity of each of these 



peaks is, however, very different: for the photoionization 
peak it amounts to 16.9 a.u.; for the first, second, and 
third ATI peaks (the latter is not shown in the figure), the 
maximum amounts to 0.31, 0.0046, and 0.00013 a.u., re- 
spectively. Clearly, the ATI peaks are weak for the short 
pulses of the considered high frequency and intensities. 

In Fig.[3l we compare the photoionization spectra com- 
puted in the present work by direct propagation of the 
electron wa ve p ackets with those published in our pre- 
vious work |10| for the same r = 3 fs pulses. One can 
see that the presently computed spectra (right panel) are 
in a very good agreement with the spectra computed by 
a totally different theoretical approach (left panel). In 
the latter approach the amplitudes of the populations 
of a restricted relevant subset of stationary states of the 
field- free system are propagated (see also above). Indeed, 
for every peak intensity indicated in the figure, each pair 
of computed photoionization spectra possesses an equal 
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FIG. 3: Photoionization spectra of the hydrogen atom ex- 
posed to Gaussian-shaped pulses of r = 3 fs duration, carrier 
frequency oi lj — 53.6057 eV, and different peak intensities 
(indicated near each spectrum) computed by two different 
theoretical approaches. Left panel: by expanding the total 
wave function in the field-free stationary states (results from 
Ref . [lOl ) • Right panel: by direct propagation of electron wave 
packet (present results, first peaks shown in Fig. [2]). 



number of oscillations caused by the dynamic interfer- 
ence, and the energy positions and relative heights of the 
multiple peak structures are very similar. 

As shown above and discussed in the preceding sec- 
tion, the spatial grid required to directly propagate wave 
packets created by the longer t = 10 fs pulses is much 
larger than for the r = 3 fs pulses. Although the prop- 
agation time required is also much longer (6r = 60 fs) 
for the longer pulse, the calculations are still feasible by 
the method and code of Ref. [3^ . In Fig. |4] we compare 
the photoionization spectra of H computed for r = 10 fs 
pulses at different peak intensities by the two different 
approaches, i.e., by the present approach (right panel) 
and by the method of our previous work [l3| (left panel) . 
The results of the latter calculations are also new and not 
available in Ref. [lO| . From Fig. |4] one can see that the 
numerically exact results of the direct wave packet prop- 
agation again agree very well with the spectra computed 
by our previous approac h [lOll . 

It is well known [38l . |39| that in strong laser fields 
each electronic state experiences an AC Stark shift, in- 
cluding the ground state, Rydberg levels, and electron 
continuous states. Since an electron in the continuum 
interacting with an oscillating field cannot have an en- 
ergy smaller than the average energy of its quiver-motion 
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FIG. 4: The same as in Fig. [3] but for r = 10 fs pulses. The 
spectra shown in the left panel are not available in Ref. [10], 
but are computed here by the method of this reference. 



(known as ponderoniotive energy), one argues that a 
strong field shifts the ionization threshold by this pon- 
deromotive energy. There is a common belief [38| that 
in the high-frequency limit, when the carrier frequency 
is much larger than the field-free ionization potential of 
a system (w ^ IP), each electronic state acquires the 
same AC Stark shift given by the universal formula for 
the ponderoniotive potential Up = Eq/Alu'^ (all quantities 
in atomic units). For this reason, at all intensities during 
the pulse the AC Stark shifts of the ground state and of 
the ionization threshold are expected to compensate each 
other [38| . In turn, this would of course eliminate the dy- 
namic interference effects as these vanish if the shift in 
the spectrum vanishes [10|. 

The results obtained and the very good agreement be- 
tween the spectra in Figs. [3] and H] computed by two very 
different approaches demonstrate that the above qualita- 
tive argumentation based on simple general assumptions 
is not correct. The present numerically exact calculations 
clearly show that the AC Stark shift of the ground state 
is not compensated by the ponderoniotive shift of the 
ionization threshold. Such a compensation is probably 
possible only for very high carrier frequencies at which 

both shifts vanish anyway, since Up > 0. It is hence 

clear now that if the parameters of the pulse, i.e., shape, 
intensity and frequency, are such that the total photo- 
electron spectrum experiences a shift, this shift results 
from the fact that these two competing effects do not 
compensate each other, and then the dynamic interfer- 
ence predicted in [T^ takes place. 

In spite of the very good agreement, a slight quantita- 
tive disagreement between the results of the two sets of 
calculations is, however, evident from Figs. [3] and 2] and 
needs to be discussed here. It is due to the incomplete- 
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FIG. 5: (Color online) Time-evolution of the cycle-averaged 
mean kinetic energy of a free electron wave packet exposed 
to a r = 3 fs pulse of frequency a; = 53.6057 eV and peak 
intensity of /o = 7 x 10^^ W/cm'^ (open circles, see text for 
details). Each circle represents the mean value obtained for 
an individual optical cycle. For comparison, the pulse inten- 
sity envelope g^{t) is also shown on the arbitrary scale (solid 
red curve). Note that the maximal value of the kinetic en- 
ergy gained during the pulse is ~ 1.2 eV at t = 0. This 
value is very different from the value of the ponderomotive 
energy Up — Sq/Alu'^ ~ 3.5 eV expected classically in the 
high-frequency limit for a field with an intensity as that at 
the pulse maximum. 



ness of the basis set of stationary field-free states used to 
expand the total wave function in [10[. This expansion 
was restricted to the ground state, Rydberg states n£, 
and photoionization ep continuum states. The conver- 
gence of the solution with respect to the included Ryd- 
berg states of different n and £ has been insured there. 
However, transitions between continuum states ei, which 
are responsible for the formation of ATI peaks in the final 
energy spectrum of the emitted electrons, were neglected 
in our previous calculations [lO|. These transitions are 
naturally included in the present numerically exact cal- 
culations. We would like to remind that even for the 
largest field intensity considered here, already the first 
ATI peak is almost two orders of magnitude weaker than 
the main photoionization peak. These weak ATI pro- 
cesses, neglected in [10|, are responsible for the slight 
disagreement between the results of the two sets of cal- 
culations. 

To better understand the impact of the pulse on the 
ejected electron, let us now investigate the motion of a 
free electron in strong laser fields. In the high-frequency 
limit [S^, it is expected that a free electron gains the 
ponderomotive energy Up, which is an average energy of 
its quiver-motion induced by a strong oscillating exter- 
nal field. In order to check this prediction in the case 
of short pulses, we start with the electronic wave packet 
given by the Is ground state of hydrogen, and propagate 
this wave packet without the attractive Coulomb poten- 
tial exerted by the nucleus. Without the pulse, this wave 
packet expands in r-space as time proceeds, but it re- 
mains unchanged in fc-space. In the presence of a pulse, 
the momentum distribution of the electron wave packet 



changes as well, and it stabilizes as the pulse expires. 

In the numerical calculation, we used a Gaussian- 
shaped pulse of r = 3 fs duration, w = 53.6057 eV 
frequency, and peak intensity Iq = 7 x 10^^ W/cm^. 
At t ~ — 3r = —9 fs, the initial wave packet is set to 
the H(ls) ground state function, which is then propa- 
gated in the presence of the pulse to i = 3t = 9 fs. 
During the whole propagation time, we have computed 
the electron kinetic energy distribution (out of the ra- 
dial momentum density distribution) at 10 time-points 
for each optical cycle, and obtain, thereby, the cycle- 
averaged energy distribution \^'^°-{e,toc)\^- This average 
distribution is then used to evaluate for each optical cy- 
cle the mean kinetic energy of the electron wave packet 
{eitoc))=Je\^^^{e,toc)\^ds. 

The mean value of the cycle-averaged kinetic energy 
obtained as described above for the free electron wave 
packet is depicted in Fig. [S] as a function of time by open 
circles (each circle corresponds to an individual optical 
cycle). Clearly, the mean kinetic energy of a free elec- 
tron increases as the pulse arrives, and it decreases again 
as the pulse expires. The final mean energy after the 
pulse is off is equal to the initial mean energy of the elec- 
tronic wave packet before the pulse was on. Moreover, 
the time evolution of the energy follows the pulse inten- 
sity envelope g^{t), which is also shown in the figure by 
a solid red curve to guide the eye. Very important, the 
maximal shift of the mean kinetic energy during the pulse 
is about ~ 1.2 eV, and it is much smaller than the value 
of Up = Eq/Auj'^ « 3.5 eV expected in the high-frequency 
limit for a field with an intensity as that at the pulse 
maximum |38| . 

This fact clearly demonstrates the failure of using the 
high-frequency approximation for a photon energy of 
~ 50 eV even for a free electron which does not expe- 
rience the attraction by the nucleus. The ponderomotive 
energy in Fig. [5] computed for a free electron wave packet 
(the field-free kinetic energy of 13.6057 eV must be sub- 
tracted from the plot) is expected to be larger than that 
for an electron which is moving in the field of the nucleus 
because of the attraction it experiences. This is only a 
qualitative argument because in the presence of the at- 
tractive Coulomb potential, the 'ponderomotive motion' 
of the photoelectron will be a part of the entire photoion- 
ization process, and its contribution to the kinetic energy 
of the emitted photoelectron will not be separable. 

As a final point, we would like to discuss which frac- 
tion of the electronic wave packet remains bound to the 
nucleus after the pulse has expired. As one can see from 
Fig. [U a part of the wave packet indeed remains bound 
to the nucleus (mainly in the ground state of the atom) . 
How to estimate the time-evolution of the population of 
the ground state? How fast is the hole created during 
ionization? In order to answer these relevant questions, 
we turn to our analytical model developed in [10| , which 
has been proven here to be reliable. As explicitly demon- 
strated there, the strong pulse induces a time-dependent 
ionization rate Tph{t), which describes the losses of the 



population of the ground state by ionization into all final 
continuum states. The ionization rate can be computed 
as the product of the total photoionization cross section 
at the chosen photon energy o-*°j*(i^) and the photon flux 
given by /(i)/w ^.E^: 



Tphit) 



'ph 



m/u. 



(9) 



The rate ^ can further be factorized as Tph (t) = Tg^ (t) . 
The rate follows the pulse intensity envelope g^{t) and 
is the strongest at the pulse maximum. The time- 
independent parameter explicitly reads F = 2Tr\duj£o/2\'^ 
[17| , where d,^ is the dipole transition matrix element for 
the ionization of the ground state computed at the chosen 
photon energy lo. 

With the help of this single parameter, one can com- 
pute the time evolution of the population of the ground 
state via the simple analytic expression [10[ 



l«/Wr=e 



2 - ^~rf'_^g'-(t')dt' 



(10) 



This analytic result was found to be in excellent agree- 
ment with the results of the full numerical calculations 
[lO| . For the Gaussian-shaped pulses considered here, 
the final population remaining in the ground state after 
the pulse has expired can be estimated analytically as 



|a/(oo)P 



.-rrj^- 



We now can try to introduce the 



time T needed for the pulse to create the hole as fol- 
lows. A pulse can ionize only a fraction of all atoms. 



which is given by iV^(oo) 



^ l^e-^I^^9'{t)dt^ 



When 



the pulse is over, this fraction defines the hole the pulse 
has created. At any time t, the population of the hole is 



Nh{t) 



1 



-r/!^3^(t')dt' 



Following the general con- 



cept of a lifetime, we define T as the time at which the 
not yet populated portion of the hole iV^(oo) — Nh{T) 
is 1/e of the final population of the hole Nh{oo), i.e., 
[Nh{oo) - Nh{T)] /Nh{oo) = 1/e. It immediately follows 
that 



T 1 - In 
g^{t)dt= 



l + (e-l)e 



-r/rL/(t)dt' 



, (11) 



from which T is easily computed. 

For strong pulses which ionize essentially the whole 
ensemble of atoms (i.e., Nfi(oo) ~ 1), condition (fTTj) 
simplifies to /_ g^{t)dt ~ l/F which is a very appealing 
result. The simple analytic expressions (|9l-|TT|) allow one 
to compute the fraction of atoms ionized during the 
pulse and to estimate the time T needed for the pulse 
to create the hole. Since the beginning of the pulse is 
difficult to define, it is convenient to consider this time 
T relative to the pulse maximum which is in our case 
at i = 0. The time T can thus be negative. For the 
T = 3 fs Gaussian pulses used in this work, for instance, 
T is equal to -H0.27, -0.19, -0.59, and -0.92 fs for the 
peak intensities 1 x 10^^ W/cm^ to 7 x 10^^ W/cm^ of 
the spectra depicted in Fig. [21 respectively. 



IV. CONCLUSIONS 

The time-dependent Schrodinger equation for the hy- 
drogen atom exposed to coherent intense high-frequency 
short laser pulses is solved numerically exactly by di- 
rectly propagating electron wave packets in space and 
time. The propagation is made without employing a 
complex absorption potential at the spatial grid bound- 
ary, and thus requires the use of very large grids even for 
the presently considered short pulses. In order to solve 
this technically challenging problem we make use of the 
efficient many-electron MCTDHF code developed in [Sj] , 
which was additionally optimized for the presently stud- 
ied single-electron problem. 

The presently computed electron energy spectra con- 
sist of a main photoionization peak and a sequence of 
above threshold ionization (ATI) peaks separated by the 
photon energy uj. Each ATI peak exhibits a pronounced 
interference pattern which resembles the multiple-peak 
structure observed in the photoionization peak which is 
due to dynamic interference |10| . For the main photoion- 
ization peak, the present numerically exact calculations 
reproduce the results of our previous calculations [lO[ 
performed by a conceptually very different theoretical ap- 
proach [iMil. 

The explicit findings allow us to conclude that the 
individual AC Stark shifts of the ground state and of 
the ionization threshold are far from compensating each 
other, as one would naively assume in the high-frequency 
strong- field limit [3^. The AC Stark effect in the elec- 
tronic continuum, as well as the dynamic interference 
effect are found to be rather pronounced. We have also 
studied the evolution of the wave packet of an electron 
which is exposed to the pulse, but does not interact with 
the nucleus. These results shed light on the quantum 
ponderomotive energies this free electron experiences in 
the short pulse. We suggest that in the real system the 
ponderomotive motion of the photoelectron in the con- 
tinuum is part of the entire photoionization process, and 
in strong short pulses its contribution to the energy of the 
emitted photoelectron is not separable from other effects. 

We would like to conclude with the following remark. 
In the course of the propagation during the pulse, the 
electron wave packet spreads over large spatial region (see 
Fig. [IJ . The parts of the total wave packet emitted on the 
rising and falling sides of the pulse are always separated 
in r-space and do not meet, as one may naively expect 
for interference to occur. Nevertheless, these two parts of 
the wave packet always overlap in A;-space giving rise to 
dynamic interference. To measure this interference, one 
should not perturb by the measurement the evolution of 
the wave packet in a rather large portion of space. 
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